Strong deflection limit of black hole gravitational lensing with arbitrary source 

distances 
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The gravitational field of supermassive black holes is able to strongly bend light rays emitted 
by nearby sources. When the deflection angle exceeds n, gravitational lensing can be analytically 
approximated by the so-called strong deflection limit. In this paper we remove the conventional 
assumption of sources very far from the black hole, considering the distance of the source as an 
additional parameter in the lensing problem to be treated exactly. We find expressions for critical 
curves, caustics and all lensing observables valid for any position of the source up to the horizon. 
After analyzing the spherically symmetric case we focus on the Kerr black hole, for which we present 
an analytical 3-dimensional description of the higher order caustic tubes. 
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I. INTRODUCTION 

Black holes have always attracted the interest of all 
physicists who hope to see General Relativity at work in 
a completely non-perturbative regime, outside any Post- 
Newtonian expansion. Since most of the information we 
receive from black holes and their surroundings is in the 
form of electromagnetic waves, one of the fundamental 
problems to be faced is the propagation of such waves 
in a black hole spacetime. The situation can be possi- 
bly complicated by an accretion disk formed by neutral 
plasma [IJ. For wavelengths at which observations are 
typically lead, the geometrical optics approximation pro- 
vides a very robust description for the propagation of the 
light rays, defined as the lines orthogonal to the wave- 
fronts. In all situations in which the plasma physics has 
little effect on the rays trajectories, the light rays sim- 
ply follow null geodesies. Then all questions that involve 
the propagation of an electromagnetic signal require in- 
tegration of the null geodesies equation. In the case of 
Kerr metric, the null geodesies have been expressed by 
Carter in terms of first integrals through the separation 
of the Hamilton- Jacobi equation Then, the integrals 
involved in these geodesies can be solved in terms of ellip- 
tic functions (see e.g. [l,|3l)- Although these analytical 
solutions are not particularly illuminating by themselves, 
they can be successfully employed to build fast and ac- 
curate numerical codes MIES, by which one can get 
particular answers, such as the shape of the iron K-line or 
the appearance of the accretion flow into the black hole 
in a future hypothetical high- resolution image 0, B Q ■ 

The problem of finding the null geodesies connecting 
a source to an observer in a curved background is usu- 
ally referred to as gravitational lensing. It has been 
pointed out by many authors that in a black hole space- 
time there are infinitely many possible trajectories for the 
photons emitted by a point-source to reach the observer 
E M, [H [H, [H, [H, For each of these trajecto- 



ries the observer will detect an independent image of the 
original source. The images can be classified according 
to the number of loops performed by the correspond- 
ing photons around the black hole. One starts from the 
primary and secondary image, which are formed by pho- 
tons performing no loops. These are already present in 
the classical weak deflection limit of gravitational lensing. 
Besides these, there are two infinite sequences of higher 
order images, formed by photons performing one or more 
loops around the black hole before reaching the observer. 
These images are progressively fainter and appear closer 
and closer to the apparent shadow cast by the black hole 
on the sky. 

The higher order images contribute much less than the 
primary and secondary images to the total flux and are 
often completely masked in situations in which they are 
not separable from the main images Therefore, the 
best chances to observe higher order images are offered by 
a black hole as massive and close to us as possible, so that 
its apparent angular size is the largest. The natural can- 
didate is the black hole at the center of our Galaxy, iden- 
tified with the radio-source Sgr A*. This is believed to 
be a 3.6 x IO^Mq supermassive black hole slowly accret- 
ing material from the surrounding environment ^is'. 
Its distance is Dql — 8kpc, so that the Schwarzschild 
radius Rsch of this black hole spans an angle of roughly 
9 /ias in our sky. A resolution of this order of magnitude 
is needed to detect the higher order images, besides the 
requirement of negligible absorption in the wavelength of 
the emitted signal by the surrounding material. In spite 
of these difficulties, the detection of higher order images 
of sources around Sgr A* should be at hand of future 
interferometers operating in the sub-mm range, where 
one expects to detect higher order images of bright spots 
on the accretion disk [9|, or in the X-ray band, where 
Low Mass X-Ray Binaries and other sources are active 

dliailil,!!!. Such images would be invaluable witnesses 

of the strong gravitational field around the supermassive 
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black hole at the center of our Galaxy; their direct ob- 
servation would thus be of striking importance in the 
confirmation of our gravitational theory. 

In addition to the considerations about their impor- 
tance, higher order images can boast two advantages 
with respect to lower order ones: if the black hole has 
non-negligible spin, they can easily form large arcs and 
additional images because their caustics have a larger an- 
gular extension compared to the first caustic; secondly, 
they enjoy a relatively much simpler analytical descrip- 
tion. 

The treatment of higher order images can take ad- 
vantage of the fact that the deflection diverges logarith- 
mically when the impact parameter reaches a minimum 
value. Then, the higher order images can be obtained by 
a simplified lens equation where the deflection terms are 
replaced by the first terms of their expansions in terms 
of the impact parameter. This procedure is analogous to 
the weak deflection limit but sets its starting point in the 
opposite regime and is thus conventionally called strong 
deflection limit. It was firstly used by Darwin in 1959 for 
the Schwarzschild black hole [IJi] and then revived sev- 
eral times 0, [H, [13, [3| until it was generalized to all 
spherically symmetric black holes ^23^ . This method was 
then applied to several interesting black hole metrics, also 
motivated by string and alternative theories [H, . 
The time dela y ca lculation for higher order images was 
done in Ref. [23|. The method was recently extended 
to the presence of external shear fields in a setup anal- 
ogous to the Chang & Refsdal lens The extension 
to Kerr metric has required several steps, starting from 
the purely equatorial case [1^ to the case of generic tra- 
jectories with equatorial observers [2lj, and finally to the 
general case |22| (in the latter two works, the treatment 
is however limited to the second order in the black hole 
spin) . An application to the Kerr-Sen dilaton-axion black 
hole has also been performed [s^. Recently Iyer and 
Fetters have found an alternative expansion parameter 
that significantly reduces the discrepancy between the 
strong deflection limit and the exact deflection formula 
[3l| . They have also explored the next to leading order 
terms in the strong deflection expansion. 

The strong deflection limit allows a simple analytic in- 
vestigation of the gravitational lensing properties of any 
black hole metric in a well-defined limit, in which the re- 
sults are easily comparable from one metric to another. 
For the Kerr metric it has been able to provide the first 
analytical formulae for the caustics and the critical curves 
involving higher order images and for lensing of sources 
near caustics. However, in the formulation used up to 
now, it has just been applied to sources very far from 
the black hole, so that the limit Dls ^ Rsch has been 
taken in all equations (here D^s is the distance of the 
source from the black hole and Rsch — 2GM/c^ is the 
Schwarzschild radius of the black hole) . 

The purpose of this paper is to remove the limitation 
to very far sources, enlarging the investigation of gravi- 
tational lensing in the strong defiection limit to sources 



placed at arbitrary distances from the black hole. We 
will thus be able to discuss the mathematical structure 
of the lensing problem and all the lensing observables 
taking D^s as an additional parameter. We will show 
that the strong defiection limit is well-defined even for 
sources inside the photon sphere, so that our discussion 
can be safely pushed up to sources lying just outside the 
horizon of the black hole. Similarly, in the Kerr metric 
we will be able to describe the caustic hypersurface from 
infinite radial distance up to the horizon. 

This paper is structured as follows. Sect. lU contains 
the new outline of the strong defiection limit for spher- 
ically symmetric black holes with arbitrary source posi- 
tion; it analyzes the lens equation and observables and 
discusses the Schwarzschild metric as a simple example. 
Sect, mil contains the extension of Kerr black hole lensing 
to arbitrary source distances, with some details moved to 
the appendix. Section [TVl contains the conclusions. 

II. SPHERICALLY SYMMETRIC BLACK 
HOLES 

In this section we shall present an updated version of 
the method outlined in Ref. [1^. Besides including the 
finiteness of source and observer distances from the black 
hole, we also make some more slight refinements that al- 
low further generalization of the method and more phys- 
ical insight. We stress the importance of the study of 
spherically symmetric black holes as propaedeutic to the 
investigation of the Kerr metric, to be tackled in the next 
section. 

As in Ref. [1^, we start from three basic assumptions 
on the spacetime metric: 

a) The spacetime is stationary and spherically sym- 
metric, so that the line element can be written in the 
form 

ds^ = A{r)dt^ - B{r)dr^ - C(r) (d-ff^ + sin^ Mc/)^) . (1) 

b) We assume that our metric is asymptotically fiat, 
so that for r — > c» we have 

A{r) 1, B{r) 1, C(r) ^ r^. (2) 

c) Furthermore, we assume that the function 
C{r)/A(r) has one minimum at > 0, corresponding 
to the radius of the photon sphere |32j . 

In some gravitational theories, the photons do not fol- 
low geodesies of the background geometry, but the self- 
interaction makes them follow geodesies with respect to 
some effective metric [25i] . These particular cases can fit 
into our treatment, provided that one uses the effective 
metric felt by the photons. Assumption (b) can be gener- 
alized to spacetimes that are asymptotically conformal to 
fiat, thanks to the conformal invariance of null geodesies. 
In this way one can include e.g. black holes with a cos- 
mological constant [1^. 

Let us give a closer look at assumption (c) and specifi- 
cally at the function C/A. We can note that asymptotic 
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C/A 



with the specific energy E and angular momentum J 




FIG. 1: The function C{r)/A{r) plotted for the Schwarzschild 
metric = 1 — Rsch/r, C = r^). In this case the static limit 
is Ts — Rsch and the photon sphere radius is rm ~ l.bRsch- 
The light-gray region is forbidden because it corresponds to 
imaginary r. The dark-gray region lies below the horizon of 
the black hole and is uninteresting for gravitational tensing. 



flatness drives the function C/A to approach at very 
large r. Moreover, in any metric admitting a static limit 
rs, such that A{rs) = 0, the function C/A diverges at 
r^. So, by continuity it must have at least one minimum 
greater than rg. As a useful reference, in Fig. [Jl we 
plot the function C/A for the Schwarschild spacetime, 
where it assumes the form /{r — Rsch)- Following Ref. 
[S^ . all stationary points ot C/A are technically photon 
spheres. However, maxima are not significant for gravi- 
tational lensing, since they are only accessible to locally 
emitted photons. We are actually interested in minima 
of C /A, since they give rise to logarithmic divergences in 
the deflection angle, as it will be clear later. Although 
it is relatively easy to find metrics which also develop 
a maximum in C /A (e.g. Reissner- Nordstrom with su- 
perextrcmal charge), it seems difficult to imagine a met- 
ric developing a second minimum. Such a metric would 
require a quite exotic source to be sustained as a viable 
solution of gravitational equations. Through assumption 
(c), we are discarding these too problematic spacetimes 
and stick to more reasonable metrics with only one min- 
imum for C/A. 

Now let us start the calculation of the photon deflec- 
tion. The spherical symmetry allows us to choose the 
equatorial plane as the plane where the entire motion of 
the photon takes place, so that 13 — ^/2 and — 0, where 
the dot denotes derivative with respect to the affine pa- 
rameter. 

The dynamics of the photon can be derived from the 
Lagrangian (see e.g. [11] for a complete discussion of null 
geodesies in Schwarzschild and Kerr spacetimes) 
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(3) 



The coordinates t and 4> are cyclic so that their conjugate 
momenta are constants of motion. They can be identified 



E 

J : 



= A{r)i 
C{r)<i>. 



(4) 
(5) 



We can choose the orientation of the polar axis so that 
J > 0. 

The last constant of motion comes from the fact that 
the photon moves along null geodesies of the metric ([T]) 
so that g^yx'^x^ = 0. From this equation, we derive the 
expression of r: 



(6) 



The angle formed by the spatial components of the 
photon momentum = (f, (p) with a normalized vec- 
tor tangent to a sphere centered on the black hole = 
(0,1/VC) is 




\p\\t\ 




(7) 



So, in any point of the photon trajectory, the knowl- 
edge of the value of the combination 



J 

E' 



(8) 



allows us to calculate the angles formed by the photon 
with respect to the radial and tangent directions. The 
photon moves radially if m = and tangentially in points 
such that u = ^ C/A. It is also easy to prove that for 
those photons reaching the asymptotic flat region, this 
quantity equals the impact parameter, defined as the dis- 
tance between the black hole and the asymptotic trajec- 
tory followed by the photon. On the basis of its immedi- 
ate connection with the observed direction of the photon, 
we will eliminate J in favor of u, where possible. 

Inversion of the radial motion can occur only at the 
points that make the argument of the square root vanish 
in Eq. ([5]) , which correspond to points of instantaneous 
tangential motion by virtue of Eq. ([7]). However, as- 
sumption (c) states that the function C/A has a single 
minimum at r^. So, a quick look at Fig. [1] convinces 
that the equation 



C(ro) _ . 
has real roots only if w > with 



(9) 



(10) 



where we have introduced the short notation Am = 
A(rm) and similarly for B and C. 

At this point it is convenient to distinguish the case 
E)hs > Tm (source outside the photon sphere) from the 
case Dls < Tm (source inside the photon sphere). We 
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shall finally find that they both lead to the same expres- 
sion for the deflection of a photon in the strong deflection 
limit, given by Eq. ([5T|) . The next two subsections deal 
with the details of the calculations in the two mentioned 
cases. 



A. Source outside the photon sphere 

Let us analyze all possibilities for the radial motion of 
a photon emitted by a source outside the photon sphere. 

Some photons will leave the source with positive r. 
Since the photons are emitted at D^s > fm, these pho- 
tons never meet inversion points and run towards the 
asymptotic region without experiencing any effective de- 
flection by the black hole. If Dql > D^s some of them 
reach the observer and give rise to the primary image, 
which is not the subject of our analysis anyway. 

Some other photons leave the source with negative r. 
If Dql < Dls, there is still the possibility that some 
of them can reach the observer without inverting their 
motion and form the primary image. But some other 
photons do not meet the observer and run towards the 
black hole. The photons with u < Um will inexorably 
sink into the black hole, since Eq. ([9]) will admit no 
real roots. If m > the photons invert their motion 
at the largest root of Eq. ([5]), which we indicate by tq 
and identify with the closest approach distance. After 
the inversion in the radial motion, these photons go back 
towards the asymptotic region and eventually reach the 
observer, giving rise to the secondary and higher order 
images. 

All these considerations can be summarized by saying 
that light rays shot at too small impact parameters are 
swallowed by the black hole, whereas those shot at larger 
impact parameters are just deflected, the limiting value 
of u being Um- Our objective is to quantify the deflection 
of these photons as a function of u. 

The azimuthal shift of the photon is 



= 





To . 


Do 






-I 








h 


Dls 


ro 



—dr, 
r 



(11) 



where we have separated the motion of the photon into 
approach phase (with r running from ZJ^g to the inver- 
sion point To) and departure phase (with r running from 
ro to Dql)- In the first integral we use the expression for 
r with the minus sign, in the second integral we use the 
expression with the plus sign. Using Eqs. ^ and ([6]), 
and reversing the extrema in the first integral, we have 



(12) 

with the short notation i = O, S and D^i — Dql, Dls 



i\(Pi — / u\ -pz \ —r — u 



Note that the integrand diverges at rp, which is defined 
as the largest root of the last factor. In order to study the 
character of the divergence, it is opportune to perform a 
detailed analysis of the function 



i?(r, u) 



C{r) 
A{r) 



(14) 



which governs the divergence of the integrand in the lower 
extremum. From the previous discussion, we know that 
R(r,u) has a minimum at for any fixed value of m; it 
vanishes at {rm,Um) by the definition of u™; it vanishes 
at (ro, u) by the definition of ro. These properties can be 
formalized by the equations 



dR 



dr 

R{rra,Uyn) = 
i?(ro,M) =0. 







(15) 

(16) 
(17) 



Since we are interested into those trajectories whose 
inversion point is very close to the minimum rm, we define 
a parameter 5 by the equation 



(18) 



Correspondingly, also the impact parameter must be very 
close to the minimum. We thus define the parameter e 
by the equation 



u„(l -I- e). 



(19) 



Then, inserting HH]) and into Eq. HT]) and ex- 
panding to the lowest order in 5 and e, we have 



= R{rm, u 
Id^R 



dR 
dr 



2 9r2 



OR, 
du 



)ume. (20) 



The first two terms vanish because of Eqs. (|T5)) and ([Tf 
and we are left with a simple relation between 6"^ and 
e, which tells us how much the inversion point ro differs 
from the photon sphere radius r^n , when we increase the 
impact parameter of the photon from the minimum value 
Um to u. Given the form of R{r, u) it simply reads 



2ul 



with 

- la^. ^ 2 



C" A 



j\2 



(21) 



(22) 



where the prime denotes derivative with respect to the 
argument and the subscript m means that the result must 
be evaluated at r^ as usual. Thus we obtain that e is of 
the same order as S^. 
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Let us analyze the behavior of R(r, u) when r is very 
close to tq. We introduce the parametrization 



1 



1] 



(23) 



with < rj < 1, and study R{r{ri), u) for small values of 
rj. Expanding to the second order in 77, S, and first order 
in e, we find 

OR 

R{ri, u) ~ R{rm,Urn) + ryn — {rm,U,n){5 + T] + 1]5 + rf) 



1 2 d'R 
--r 

2 " 



dr 



(24) 



Again, the first two terms vanish because of Eqs. HSj) 
and Moreover, the e-term cancels with the remain- 
ing (J-^-term because of Eq. ((20)) . The behavior of i?(77, m) 
for small rj is thus 



i?(77,u) = (2577 + 772) + 0(7^2). 



(25) 



Now that we have found the dominant terms in R{ri, u) 
for 77 close to 0, which corresponds to r close to rg, we 
can return to the integral in Eq. ([13]). Changing the 
integration variable from r to 77, the integration ranges 
become [0,77^], where 77^ = 1 — r^/DLi. Each integral 
assumes the form 




[Riv.u)]-^'^ 



(1-77)^ 



rdrj. (26) 



Now we add and subtract a term containing the diver- 
gence of the integrand at small values of 77. We then 
separate the integral into two parts: 



(27) 
(28) 




C(77)™ (1-77)2 



Bn 



Pm V Cm \j25ri + 77^ 



drj. 



(29) 



Of course, the sum of the two integrals Id and In is 
just the original integral ((26l) . but now the first integral 
is elementary and reads 



In = r„ 



■log 



(30) 



The divergence for (5 — > appears explicitly in Id, while 
In is the integral of a regular function and does not di- 
verge any more for (5 — > 0. Fig. [5] illustrates an interesting 
comparison between the original integrand of Eq. 
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FIG. 2: The plot illustrates the integrand of A(f)i in Eq. ((26)) 
as a function of ry for S = 0.025 (e = 0.001) in the case of 
Schwarzschild black hole; the dashed line is the integrand of 
lo in Eq. (|28|l . the difference between the two curves being 
the integrand of In in Eq. . 



and the integrand of lu in Eq. ([28)) taking the special 
case of Schwarzschild metric as an example. Although 
the integrand of lu is a drastically simplified form of the 
original one which approximates it for 77 very close to 
zero, it turns out to be a very good approximation in the 
whole range of 77. Only for 77 close to 1 we see a sensi- 
ble difference. Such difference is stored in the integrand 
of /fl. As b tends to zero, the contribution by Ju be- 
comes more and more dominant with respect to that of 
Ib. as the divergence of the integrand becomes stiffer and 
stiffer. These considerations are a good premise to the 
strong deflection limit. 

Until now we have done no approximation. We have 
just added and subtracted some terms and made some 
changes of variables. The expression ([27]) is still exact. 
The strong deflection limit amounts to save the first dom- 
inant terms in the expressions for Id and /ij as (5 0. 
We have a logarithmic divergent term in Id and then 
some terms converging to constant values. We also note 
that the parametrization (j23p tends to 



(31) 



1-77' 



when (5 — !■ 0. Consequently, the integration limits 77^ tend 
to 77, = 1 - r,n/DLi. 

After the truncation of the expansion in (5, we have 



A(^, =alog^ + 6, + o(<50), 




(32) 
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where the coefficients a and bi are given by 



B„ 



(33) 



bi = J gi{r])drj ^ J Sign{r]) gi{r])drj (34) 





C(r;) 



(1 - r,)2 




(35) 



where gi{i]) is just the integrand of Ir in the Hmit for 
(5^0, also implying u Um- The function Sign{ri) has 
been introduced only for uniforming the expression of bi 
to the corresponding terms we shall derive in section fll BI 
in the case of a source inside the photon sphere. For a 
source outside the photon sphere, the variable rj is always 
positive. Finally, the deflection suffered by the photon is 
quantified by the full azimuthal shift A(j), given by 



A(j) — a log 



6^ 



bs- 



(36) 



When the source and the observer are very far from 
the black hole, it makes sense to define a defiection angle 
as the difference between the azimuthal shift suffered by 
the photon minus tt, which represents the total azimuthal 
shift of a photon travelling in a flat space without the 
black hole on a rectilinear trajectory. This concept be- 
comes ill-defined for sources and observers that are not 
in the asymptotic flat region. 

The fact that source and observer are at finite distances 
is encoded in the presence of 770 and 775. Setting them 
to 1, the defiection angle so derived coincides with the 
expression originally given in Ref. [2^ . 

With the definition of (3m given by Eq. (22) and Eq. 
(I14p . we can formulate an explicit expression for the co- 
efficient of the logarithmic term in terms of the metric 
functions 



n" A — C A" 



(37) 



This coefficient is independent of the source and ob- 
server positions. 



B. Source inside the photon sphere 

Now let us consider the case in which the source is 
inside the photon sphere, but still outside the horizon. 
The photons leaving the source with negative f sink into 
the black hole. As for the photons leaving with positive 
f, we have two possibilities: those starting with u > 
Um meet an inversion point before reaching the photon 
sphere radius rm- Therefore, they fall back into the black 



hole. The photons with u < Um rneet no inversion point 
and escape towards the asymptotic region. The observer 
will thus see the deformed images of a source inside the 
photon sphere. 

Recalling the relation between u and the angle formed 
by the photon momentum with the tangential direction 
(Eq. (O), we can re- interpret this discussion noting that 
only photons shot along the radial direction {u — 0) or 
very close to the radial direction (u < Um) will be able 
to escape to infinity. Photons shot at larger angles with 
respect to the radial direction invert their motion before 
crossing the photon sphere. It is interesting to note that 
the angle formed with the tangent direction by photons 
emitted at Dls < fm decreases until they cross the pho- 
ton sphere. After that moment, they align more and 
more with the radial direction as they move farther and 
farther from the black hole. 

At first sight, one may think that this situation is very 
different from the one described in the previous subsec- 
tion. Actually, even in this case it is possible to define a 
strong deflection limit, corresponding to photons with u 
just slightly smaller than Um- Let us see this in detail. 

The azimuthal shift of the photon is 

<t>f DoL . 

A(l) = J d(f>^ J ^dr 



DoL 



uJ-[R{r,u)]-^^^dr (38) 
V o 



Dr 



with Dls < frn < Dql- Even if the function R{r,u) 
never vanishes, it becomes minimum at r — r,,,,. Corre- 
spondingly, the integrand has a maximum at this point 
and is largely dominated by this peak at if u is very 
close to Um- So it is convenient to revisit the analysis of 
the function R(r,u). 

Now we have to keep in mind that u < Um, so that the 
parametrization ((TH]) yields e < 0. As pointed before, the 
function R(r,u) has no real roots when u < Um, and in 
fact, Eq. (PT|) with e < gives an imaginary value for 5, 
so that the inversion point tq = ^^(l -I- S) is no longer 
a real number. Moreover, it is again convenient to intro- 
duce the parametrization (j31[) , but this time extended to 
r < Tm corresponding to < 0. Thus now the 77 variable 
ranges in the interval 1 — ^ < 7/ < 1. 

We can now study the function R{r{r]),u) for small 
values of rj as in the previous subsection and find 

OR 

R{ri,u) ~ R{rm,Um) + -7^{rm,Um)rm{r] + V^) 
or 

1 d^R. \ 2 2 dR , , . . 

^2~dr^ "^'^ "^'^ - -^(^m,u„)wm|e|. (39) 



As usual, the first two terms vanish because of Eqs. (fu 
and ITBl) . As for the e-term, we can replace it by the 
corresponding 5^-term through Eq. (PO)) . We just have 
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to keep in mind that now 5"^ < 0. The behavior of R{ri, u) 
for smaU 77 is thus 



R{t],u)^(3^{~S^ + tj^) + o{7^''). 



(40) 



Returning to the integral ((38|) . we can change the in- 
tegration variable from r to ry using Eq. (|3ip to get 



A(p = 



rdrj, (41) 



where the integration extrema are already in the form 
r]i — 1 — Tml Dn. Note that, since the source is inside 
the photon sphere, we have rjs < 0. 

As before, we add and subtract a term that contains 
the main structure of the integrand, that is the peak at 

= 0, corresponding to r = r^. We have 




/3m V Cm \/ — S^ + 77' 



-.d^ 



(42) 
(43) 



^ V Cm ^-(52 + 7^2 

The first integral is again elementary and reads 



(44) 



■ log , .r — — r ' (45) 



^m/3m v^-52 + nl - |77s| ' 



where we have made the sign of 775 explicit. The second 
integral contains an integrand that is regular everywhere 
for (5^0. Fig. [3] illustrates a comparison between the 
integrand of Eq. (|4T|) and the integrand of Id in the 
Schwarzschild case. We see that even if there is no di- 
vergence for 77 = 0, the integrand has a very pronounced 
peak that dominates the integral. The peak structure 
is catched by Id, while the wings are corrected by the 
contribution oi Ir. As e — > 0, the peak grows larger and 
larger, dominating the wings. 

Now we make the strong defiection limit approxima- 
tion, by requiring that e and consequently 5"^ is small. 
Saving the logarithmically divergent term and the con- 
stant terms, we have 



— a log 



<52 



+ 5so + o(5°), 



(46) 



where the coefficient a is still given by Eq. ([33|) and bso 
is 



vo 

bso = J 9iiv)dv, 



(47) 




-0 . 4-0 .2 0.20.40.60.8 1 

n 

FIG. 3: The integrand of A<f> in Eq. (|4T|l as a function of rj 
for S = 0.025i (e = —0.001), compared to the integrand of Id 
in Eq. (|43p in dashed style. 



with 31(77) still given by Eq. ([55)1 . It is interesting to 
note that the argument of the logarithm remains positive, 
since both and 775 are negative. 

We can split the integral in Eq. (|47p in two parts 



bso ^bs + bo 







ns 




J Sign{ni) giir])dr] 


ns 





■no 


•no 


bo ^ J 9iiv)dv = 


J Sign{rj) gi{rj)dr] 









(48) 



We observe that bo and bs have the same formal expres- 
sion as bo and 65. 



C. Deflection and higher order images formulas 

As a result of the previous two subsections, we have 
a unique expression for the photon deflection, which can 
be conveniently stated in terms of e, which represents the 
impact parameter shift from the minimum value (see Eq. 
pip. By Eq. §^ and 1^ we finally get 



A(j) ~ —a log ■ 



where we have defined the coefficient 
2P„ 



+ b + TT, 



b ~ a log ■ 



+ bo + bs — TT. 



(51) 



(52) 



For quick reference, e is defined in Eq. (|19p . Um is 
given by Eq. pU]) . /?m is given by Eq. a is given 

by Eq. ([33l), 77^ = 1 - rm/DLf, bo and bs are given 
by the integrals ([M)) . where 31(77) is found in Eq. ([55)1 
and R{r,u) in Eq. (fT4|l . having changed the integration 
variable from r to 77 through the parametrization pip . 
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This expression for the total deflection of the photon is 
vahd for any position of observer and source. Even if we 
have not exphcitly considered it, time-reversal symmetry 
warrants that Eq. (jSip is valid even in the unrealistic 
case of an observer inside the photon sphere (provided 
that one correctly relates e to the sky coordinates of an 
observer in a strongly curved region). The only approxi- 
mation performed is that the impact parameter u is very 
close to the critical value Um- 

The general lens equation for spherically symmetric 
black holes can be simply stated as 



in fm/Dn. Recalling that rji = \ — rm/D^i, we have 



bs — mod 27r. 



(53) 



Fixing the origin of the azimuthal coordinate in such a 
way that (j)o = tt and using the expression of the total 
deflection in the strong deflection limit ([5T|l . we can eas- 
ily solve the lens equation and find the position of the 
images. In general we have 



En = 110 Vse 



i>+<f>s- 



(54) 



where n denotes the number of loops done by the photons 
before reaching the observer and 0s € [— tt, tt]. Of course, 
the strong deflection limit becomes exact in the limit n — s- 
oo but is typically a very good approximation already for 
n = 1, as will be shown in the next subsection. 

For an observer in the asymptotic region, which is the 
most physically interesting case, the angular separation 
between the direction of arrival of the photon and the 
direction of the black hole is simply 9 — u/Dql- So, we 
have 



OL 



(55) 
(56) 



6m is usually called the angular radius of the shadow 
of the black hole, since all images of sources outside the 
photon sphere reach the observer from angles 9 > 9m and 
the region within the angular radius 9m appears empty. 
However, when the source is inside the photon sphere, 
Dls < fm and then rjs < 0. We thus have e < and 
the sequence of images will appear within the shadow of 
the black hole, with the lowest order ones closer to the 
center and the higher order ones closer and closer to the 
shadow border. 

The study of the Jacobian of the lens equation confirms 
that the critical curves are simply given by ([55]) and ([54]) 
with (^5 = for standard lensing (ps = for retrolensing. 
They are circles outside the shadow for sources outside 
the photon sphere and inside the shadow for sources in- 
side the photon sphere. The caustics are always pointlike 
and are located behind and in front of the source. Alto- 
gether, the caustics cover the whole optical axis as Dls 
varies from -|-c» to the static limit Ts- 

It is interesting to take the limit of the total azimuthal 
shift A0 for Dol,Dls ^ Tm and calculate the first order 



A(j) = —a log 



2f3„ 



DoL Dls 



gi'n)dr] 



+g{l)ivo - 1) + gWivs - 1) + oir,, - 1). 



(57) 



g{l) can be calculated using the asymptotic limit of all 
metric functions evaluated in rj that appear in its expres- 
sion and taking the limit for 77 — > 1. It is simply 



5(1) 



Summing up, we get 

A(/) — TT = a — 9m — 



(58) 



(59) 



where a is the deflection angle calculated on the asymp- 
totic trajectories [Dls = Dql = 00), Om is defined in 
Eq. dSi) and 



(60) 



Dls' 



is the angular size of the shadow of the black hole as 
measured by a distant source. The lens equation then 
becomes 

-(j>S^a~9m~9m (61) 

The first correction to the lens equation with far source 
and observer is thus universal and simply takes into ac- 
count the geometry of the lensing problem (compare with 
the discussion of the lens equation in Ref. [s^]). 



D. Testing the formulas in the Schwarzschild case 

In this subsection, we shall specify all our general for- 
mulae for black hole gravitational lensing in the case of 
the simplest possible metric. This will allow us to under- 
stand the sense, the validity and the power of the strong 
deflection approximation throughout the range of possi- 
ble source positions. 

The Schwarzschild metric reads (with G = c = 1) 



2M 



A(r) = 1 - 

r 

B{r) = (1 - — j 
C{r) =r2. 



(62) 

(63) 
(64) 



The minimum of the function C /A is at = 3M. 
Correspondingly the minimum impact parameter Um = 
\/Cm/Am = 3-\/3M as well-known [l^, . 

Now let us calculate the coefficients of the deflection 
formula in the strong deflection limit. As already noted 
after Eq. (|37p . the coefficient of the logarithmic term 
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is independent of the source and observer positions. So, 
Eq. (|37p simply gives the already known value 



a = 1. 



(65) 



This can be expected, since the logarithmic term is 
a characteristic of geodesies winding around the photon 
sphere which does not depend on the start and arrival 
point. 

The constant coefficient in the deflection formula is 



6 = -TT + 5 log[6] + bo + bs 



h = -21og 



18Af 
D 



Li 



Putting everything together, we have 
A0 — TT = — log ■ 



(66) 
(67) 

-6, (68) 



{l-MIlDLs){l-iM/DoL) 
which reduces to the well-known formula [lo| 

a = -loge + log[216(2- \/3)2] - TT (69) 

in the limit Dol,Dls ^ oo. Eq. (j68|) is the generaliza- 
tion of Darwin's formula (j69p to sources and observers at 
finite distance from the black hole, the only approxima- 
tion remaining |e| <C 1 (see Appendix \X\ for a discussion 
of formulae expressed in terms of alternative perturbative 
parameters). In order to test our formula for the deflec- 
tion of a photon, we can use it to calculate the radius of 
the critical curves. 

The angular radius of the critical rings is given by Eq. 
(151)1 with 05 = or TT. Explicitly, for the Schwarzschild 
metric we have 



DoL 



(1 + efc) 



_ 216(2~V3)(2i^LS-3) ^_,, 



(70) 
(71) 



where k is an even number in the retro-lensing case 
(05 = 7r), and an odd number in the standard lensing 
case (05 = 0). The first critical curve for fc = 1 is cre- 
ated by photons experiencing weak deflection and is be- 
yond the range of validity of Eq. ([70)1 . The critical curve 
with fc = 2 is the first retro-lensing ring, while for k — 3 
we obtain the first higher order Einstein ring of standard 
lensing. 

The displacement of the first relativistic Einstein ring 
from the black hole shadow is shown in Fig. 2^, where 
it can be appreciated how 63 stays small throughout the 
range of source distances, validating the strong defiection 
limit approximation. For D^s 3> M, the ring radius 
tends to its asymptotic value. Our analytic formula for 
the critical curve nicely joins the region within the photon 
sphere Dls < 3M to that outside the photon sphere 
Dls > 3M. The divergent term in the deflection formula 
at Dls = 3M forces e = 0, in order to keep A(f) finite. 




1000 



1x10 




FIG. 4: (a) Angular radius of the first relativistic Einstein 
ring relative to the black hole shadow for a source behind 
the black hole as a function of the source distance measured 
in Schwarzschild radii (es = {63 - 9,n)/0m, 2M = 1). (b) 
Difference between the value of 63 calculated by the strong 
deflection limit formula and the one calculated by the exact 
formula. 



Since this is true for any order k, we can conclude that 
all higher order images of a source right at the photon 
sphere collapse into one degenerate image appearing right 
at the shadow border. As Dls < 3M, e becomes negative 
as expected. This means that the sequence of Einstein 
rings for a source inside the photon sphere is reversed: 
the brightest rings will appear closer to the center of the 
shadow and the fainter will be farther, approaching the 
shadow border as k grows more and more. 

The Schwarzschild metric is simple enough to allow an 
exact integration of the azimuthal motion in terms of el- 
liptic integrals. It is thus very important and instructive 
to compare the critical curve radius calculated by the for- 
mula (|5ip obtained in the strong deflection limit to the 
exact Einstein ring position, calculated using the exact 
deflection angle. The difference between the two posi- 
tions is plotted in Fig. [41d, where it can be appreciated 
that it stays of the order throughout the range of Dls, 
testifying the accuracy of the strong deflection limit as a 
powerful approximation to describe higher order images. 

The updated analysis of extreme gravitational lensing 
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by Schwarzschild black holes presented here can be re- 
peated for any kind of spherically symmetric black holes, 
using the formulae derived in this section. It is interest- 
ing to consider that some metrics could give the same 
lensing observables for sources at infinity, whereas they 
could be distinguished in gravitational lensing of sources 
at small distances from the black hole. Therefore, the in- 
troduction of 15^5 as a new parameter enriches the arena 
for the comparison of different metrics. 



the equations 




with 



III. KERR BLACK HOLE 

In this section, we shall describe gravitational lensing 
of sources placed at arbitrary distances from a spinning 
black hole. With respect to the analysis of Ref. [22], the 
finiteness of Dql and Dls only intervenes in the calcu- 
lation of the radial integrals. As a consequence, the first 
three sections of Ref. [22], concerning the description of 
the unstable circular orbit and the shadow of the Kerr 
black hole, remain unaffected. Their content is briefly 
reviewed in the following subsection. The modifications 
in the radial integrals are reflected in the lens equation 
and its Jacobian. Consequently, also the critical curves 
and the caustics contain a dependence on the source po- 
sition. They are described in Section [ill Bl Finally, lens- 
ing of sources near caustics is updated in Section UlI CI 
Throughout this section, we shall preserve the spirit of 
Ref. [24] expanding all quantities to the second order in 
the black hole spin a. The perturbative expansion is of 
key-importance to keep all calculations fully analytic up 
to the final results. At the same time, the second order 
approximation proves to be very reliable up to values of 
the black hole spin a ~ 0.1, as noted in the comparison 
with numerical results [2 ill. 



e = Q + a^ cos^ d~J^ cot^ i9 (77) 
R= + (a^ - J2 - qy + (Q + ( J - af)r 

^a^Q. (78) 

J is the component of the angular momentum of the pho- 
ton along the spin axis and Q is the Carter integral 0], 
related to the total angular momentum of the photon (we 
set the specific energy i? to 1 by a suitable choice of the 
affine parameter). 

An observer in the position (Dql, ^Oi (t'o)i defines an- 
gular coordinates {0i,O2) in the sky, such that the black 
hole is in (0,0) with its spin projected along 92- A pho- 
ton travelling on a geodesic characterized by constants of 
motion J and Q hits a distant observer from the direction 



(79) 
(80) 



= ±D 



OL 




J2 



where /io = cos -do as usual [32 



Shadow of the Kerr blaek hole 



A. Derivation of the lens equation 



The Kerr metric in Boyer-Lindquist coordinates [36j is 



ds 



2 „ A-g^sin d _ p 



y+a^) -a^Asin^i? 



P 



2ar sin^ 

H ^ dtdcp 

P 

r — r + a , 
+ cos^ d. 



■ sin^ dd(j)^ 



(72) 

(73) 
(74) 



Distances are measured in units of the Schwarzschild ra- 
dius {2MG/c^ = 1) and o is the specific angular mo- 
mentum of the black hole, running from (Schwarzschild 
black hole) to 1/2 (extremal Kerr black hole) in our units. 
The Kerr geodesies are described in integral form by 



Among all photon trajectories ending at the observer, 
there is a family of trajectories that approach an unsta- 
ble circular orbit around the black hole when traced back 
asymptotically in the past. This family can be parame- 
terized by the parameter ^ ranging from —1 to 1. The 
constants of motion identifying the geodesies of this fam- 
ily are then given by 



3V3 



S,^l-pl-a{l-pl){l + e) 

2e-2pi{i-e)]m 



Qm(0 



3V3 



[i-{i-ph)e 



27 



-iVia^^J I - pI[1 ^ 
+3(1-M?,)Y]. 



2 

Po 



(1 - pDe] 
pDe 



(82) 



11 



The radius of the unstable circular orbit asymptoti- 
cally approached in the past is 



3 2 

— 77 



2 V3 



27%/3 



a3^(5 + 6Aio)\/l-Mo + 0(a*). (83) 



2. Strongly deflected photons 

We now introduce the following parametrization of the 
observer sky by 



(91) 



Correspondingly, a distant observer sees such photons 
from the directions 



DoLdl^r. 



,^^[5-2^,l-2e{l-t^l)], (84) 



DolO 



2,m 



Mo 



Ta2^^^[l + 2m^ - 2e{l - M^)]. (85) 



3\/3 



Eqs. (|M|) and (I85p define a curve in the observer sky 
as ^ varies between —1 and 1. This curve represents the 
border of the shadow of the black hole, in the sense that 
all photons emitted by a source outside the unstable cir- 
cular orbits reach the observer from directions outside 
this border. In the Schwarzschild limit, the shadow bor- 
der is simply a circle of radius 3\/3/2Dol, whereas for 
generic values of the spin, the shadow satisfies the ellipse 
equation 



Af 



fl2 



(86) 



One half of the sky is covered as ^ varies from —1 to 1 
and e varies from —1 to +oo. The double sign in 02, m 
selects which half of the sky we are covering. We are 
interested into strongly deflected photons, which corre- 
spond to very small values of e. One may regard e as the 
relative displacement of the photon direction from the 
shadow border. As e — > 0, the photon spends more and 
more time close to the unstable circular orbit, and per- 
forms more and more loops around the black hole before 
emerging. 

Using Eqs. (|7^ and ([50)1 we can calculate the values of 
J and Q identifying the geodesies of such photons. For 
each value of J and Q, we can calculate the inversion 
point tq in the radial motion examining the roots of the 
function i?,(r, J, Q), defined by Eq. ([75)1 . It is immediate 
to see that R{r, J, Q) satisfies the properties 



dR 
dr 







i?(ro,J,Q) = 0. 



(92) 

(93) 
(94) 



Analogously to the spherically symmetric case, we can 
define the parameter 5 by 



= r™(l + 5). 



(95) 



with the origin shifted rightward by 



2a 

DoL ' 



and semiaxes given by 

3\/3 



2DoL 
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3V3 /3V3 2yl 



(87) 

(88) 
(89) 



5 is thus the relative displacement of the inversion point 
from the unstable circular orbit radius. Inserting the 
expression of J and Q as functions of ^ and e and solving 
Eq. ((94| for 5 to the lowest order in e we find that the 
two parameters are related by 



where we have introduced the compact notation 



, (96) 



(97) 



with ellipticity 

e^l-^ = ^a2(l-M^) (90) 

For more details on the derivation of the parameteri- 
zation ([55]) and of the shadow border, see the deep and 
detailed discussion in Ref. ^22i] . 



So, for any photons hitting the observer from a direc- 
tion close to the shadow border, we can immediately find 
the corresponding value of the inversion point in its mo- 
tion around the black hole. If the photon comes from the 
interior of the shadow (e < 0), then the inversion point 
becomes complex, signaling the fact that the photon fol- 
lows a geodesic without inversion points in the radial 
motion. 
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3. From the geodesies equation to the lens equation 

Now we can come to the resolution of the radial inte- 
grals appearing in the geodesies equations ([75)1 and (ffS)) 
We indicate them by 

dr 

71 



Ii=± 



(98) 



± 



a J 



-dr. 



(99) 



These integrals must be performed along the photon 
path from the source to the observer. The double sign 
indicates that one has to sum the contributions with dif- 
ferent signs of f coherently. The task is simplified by 
the observation that the function i?(r, J, Q) determining 
the radial motion in the Kerr metric satisfies the same 
properties as the function R{r, u) determining the radial 
motion in spherically symmetric metrics (compare Eqs. 
(HHl-dMl) to Eqs. dTCJ-dUD). So, we can repeat exactly 
the same steps of the analysis of the function R{r, u), with 
the trivial extension to the presence of two constants of 
motion. The integrals are then easily solved to 

h= -ai\ogS + bi+ci{DLs)+ci{DoL) (100) 

/2= -a2\0gS + b2+C2{DLs)+C2{D0L), (101) 

where the coefficients ai, 6i, 02, and 62 remain the same 
as those given in the appendix of Ref. (2^ , whereas the 
explicit expressions of the new functions ci and C2 are 
shown in the Appendix [B] at the end of this work (Eqs. 
(|B3I) and (jB4[) ). We have stored the whole dependence on 
the source and observer distance in these two coefficients. 
They are defined in such a way that they vanish when 
their argument goes to infinity, so that the expressions 
for distant sources and observers are recovered. Although 
we have provided the expressions of the integrals for ar- 
bitrary observer distances, we shall consider Dql S> 1 
from now on for simplicity. If one is interested to study 
gravitational Icnsing with observers close to the black 
hole, it is easy to recover the relevant formulae, since the 
dependence on D^s and Dql is interchangeable, thanks 
to the symmetry Dls ^ Dql in Eqs. ^(TUU\i and (fTUT|) . 
In the angular integrals 



Ji -± 



J2 = ± 



1 

9 

CSC 



d^ 



-dd 



(102) 
(103) 



there is no dependence on the source and observer dis- 
tance. Therefore, we can safely exploit the expressions 
in the appendix of Ref. [2^ without any more concern. 

Finally, it is convenient to replace the parameter 5 
(which, as we recall, represents the relative displacement 
of the inversion point from the radius of the unstable 
circular orbit) by a new variable 



-0 = — 2 log 5 + log 



144(7-4V3)(2i?L5-3) 
(TSAE^ + ^/i + DLsY 



(104) 



From the physical point of view, ip simply represents 
the equivalent azimuthal shift of a photon deflected by a 
Schwarzschild black hole. 

Once all integrals are solved, we can rearrange the in- 
tegrated geodesies equations in the typical form of a lens 
equation, moving the source angular coordinates on the 
left hand side and leaving everything else on the right 
hand side 



A^s = lJ-s{i^, C; /-to, Dol,Dls) 



(105) 



The lens mapping relates the source coordinates {fJ.s, 4's) 
to the set of intermediate variables {tp,S,), which charac- 
terize the photon geodesic by the amount of deflection 
and its orientation in space respectively. The sky coor- 
dinates (01,02) are related to (V',^) by Eqs. ^ 
and (|104p . The observer position and the source distance 
play the role of parameters of the lens mapping. 

In the following subsections we shall put in evidence 
the main features of the lens mapping through the analy- 
sis of its critical points and the discussion of gravitational 
lensing near the critical points. 



B. Critical curves and caustics 

In order to find the critical points of the lens equation, 
one has to calculate the Jacobian determinant 



The critical points of the lens mapping are the solu- 
tions of the equation det J{ip, ^) = 0. The critical curves 
are the corresponding points in the observer sky (0i, 62) 
through Eq. ([91]l and the caustics are the images of these 
points in the (/is , (f)s ) space through the lens mapping. 
The full procedure is straightforward and is detailed in 
Ref. [l^l- Here we just state the updated results with 
the encompassment of the finiteness of I?ls- 



1. Critical points in the (V'lC) space 

As mentioned before, represents the equivalent az- 
imuthal shift of a photon moving in the Schwarzschild 
metric obtained turning the black hole spin off. It is thus 
not surprising that the zero-order critical points are sim- 
ply given by ^p^"^ ~ kir, with k being an integer number. 
k is the number of half-loops performed by a photon in 
the Schwarzschild metric. We will often refer to this in- 
teger number as the critical curve order or caustic order. 
k = 1 gives the azimuthal shift of a photon emitted by 
a source behind the black hole and weakly deflected by 
the black hole. Such photons are not the subject of our 
analysis. The first interesting case is k = 2 correspond- 
ing to photons emitted by a source in front of the black 
hole and backscattered by the black hole (retro- lensing) . 
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Photons with fc = 3 are again emitted by a source behind 
the black hole, but now the photons perform a complete 
loop around the black hole before reaching the observer. 
Summing up, odd critical orders are involved in gravita- 
tional lensing of a source behind the black hole, whereas 
even critical orders are involved in rctro-lcnsing. For each 
k we have a different critical curve, which physically cor- 
responds to the degenerate image of a source placed on 
the corresponding caustic. Starting by the zero-order so- 
lution, the first and second order correction can be found 
solving the Jacobian equation order by order. For each 
k, we have 

^k7: + aV^^O + (107) 
with ip^j}^ and ip'j^'^ given by 



4'' 



LS 



v-r =- 



{2Dls - i)VS + DLs 

^ ^(3-2^*^,-3^") + 32(1-^2^ 



(108) 



18 



3V3{2Dls - 3)2/01^(3 + Dls] 



3/2 



{2DLs-3)i3 + DLs) 



2[Vd^~VS + Dls 
•(2i?i5-3)(3-2/4_3f ) 

+26Dls - 27 - (ADls - 18)mo " (18^l5 - 27)^"^ 
-108£»Lsf I 



(109) 



In the limit 13^5 —> oo the first order correction vanishes 
while the second order correction reduces to the first row. 
These correctly reproduce the results of Ref. [11] . As ^ 
varies between —1 and 1, we recover all critical points in 
the (■0,0 space. 



which reduces to the shadow shift in the limit D^s ~^ oo. 
For generic values of the source distance, the center of 
the critical curve does not coincide with the center of 
the shadow. However, the displacement is very small, 
since it is proportional to (given in Eq. ([71]), which 
is very small in the strong deflection limit. We can also 
note that the degeneracy between the black hole spin 
a and the observer position y^l ~ fi^ that was pointed 
out in Ref. 22^ holds even when the source is at finite 
distance. Fig. [5^ shows the dependence of the shift with 
Dls. For any value of the critical order k, the shift tends 
to the shadow shift for large values of Dls- The first 
retrolensing critical curve {k — 2) is the most displaced 
one. 

The major semiaxis of the critical curve is oriented 
along the projection of the spin on the observer sky. Its 
explicit expression is 



A2.k = 



3V3 



a 



Afil - 9ck{2fil - 3) 



2DoL ' DolV 
ek \2^/lhrs^/i + DLs{'iDLs- 



3) 



4V3 

-1 



2{2Dls - 3)(/dI^ - V3 + ^Ls)'(3 - 2^ll) 
f9(2i^iS-3) + 2^2^(2I?is + 9)]}, (112) 



and is plotted in Fig. [S)d as a function of the source 
distance. We can appreciate that all critical curves are 
external to the shadow as Dls > 3/2 and internal to it 
when Dls < 3/2, as in the spherically symmetric case. 
The larger the order of the critical curve, the closer the 
curve is to the shadow, as ^ for k oo. 

Rather than giving the expression of the minor semi- 
axis, it is instructive to examine the ellipticity of the 
critical curve, defined by 



2. Critical curves 



ek 



1 



^2,fc' 



(113) 



Starting from the solution of the Jacobian determinant 
equation in the ("0,0 space, we can construct both the 
critical curves and the caustics. As for the critical curves, 
it is sufficient to use Eq. |9T|) with e expressed by Eqs. 
([M)) and (|104p in terms of tp and put = 0^ as given 
by Eq. (|107p . In this way one gets the critical curves in 
the form (0i,/c(C), ^2,fc(?))- These expressions are lengthy 
and are not very transparent. However, it is straightfor- 
ward to prove that they satisfy the ellipse equation at 
the second order in a 



71, fc 



9o,fc)^ 



42 



^2,fc 



l + o{a^). 



(110) 



The center of the critical curve is shifted by the quan- 
tity 



DoL0o,k = 2a^l - 1 



i\/iDLsek 



{2Dls - 3)V3 + Dls 



(111) 



Of course, Ai ^ can be easily derived by this expression 
if one is interested in it. To the second order in a and to 
the first order in Ck, the ellipticity is 

2 2 J 2 (27cfe - 8)efc 
Cfc =a (1 - /io) I - + - 
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4e 



V3 



■ Cfc 



9/dI^(3 + DLsf^^{2DLs - 3) 
V3{729 - 22bDLs ~ 348i:'2 
72dYs^2, + Dls]]. 



LS ^92Dls + ^8Dls) 
(114) 



This quantity reduces to the first row when Dls 00, 
which is the same as that given in Ref. [l^l safe for the 
fact that we stop here at the first order in Ck in order to be 
consistent with the strong defiection limit approximation. 
For large values of fc, efc ^ and the ellipticity of the crit- 
ical curves tend to the ellipticity of the shadow ([50]) . Fig. 
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FIG. 5: (a) Apparent shift of the center of the critical curves 
with respect to the black hole position as a function of the 
source distance; from the bottom up, the curves are for 
k = 2, 3, 4, respectively, (b) Major semiaxis of the critical 
curves with k = 2,3,4 (from top to bottom) as a function of 
the source distance, (c) Ellipticity of the critical curves as a 
function of the source distance; from the top down the curves 
are done for k = 2,3, 4, respectively. All plots are done for 
a — 0.1 and equatorial observer no = 0. 



[St shows the ellipticity of the first three relativistic crit- 
ical curves as functions of the source distance. We note 
that all curves are more elliptical than the shadow when 
Dls is large. But for sources just slightly farther than 
the unstable circular photon orbit, the ellipticity of the 
critical curves equals the ellipticity of the shadow. This 
happens at a value of D^s slightly greater than 3/2 and 
different for all critical curves. Below this distance, the 



ellipticity of the critical curves becomes smaller than that 
of the shadow. Finally, we can note that the ellipticity 
remains a function of ai/l — fj,Q, thus preserving the de- 
generacy between the black hole spin and its orientation 
relative to the line of sight. 



3. Caustics 

The caustics are obtained evaluating the lens mapping 
in the critical points ipk given by Eq. (jl07p . To the zero 
order, the Schwarzschild caustics are recovered. They are 
pointlike and placed behind the black hole for k odd and 
in front of the black hole for k even. As the black hole 
spin is turned on, the caustics drift from the optical axis 
and acquire a finite extension. Their explicit expression 



Ms 



= (-l)Vo ± RkJl - - ef^\ (115) 



0S = (1 - fc)7r - A0fe 



Rk 



with the drift given by 



Vl-Mo" 

3(2 - Vsf 

9{Dls-1) 
and the semi-amplitude given by 

Rk = a^(l - fip) [cfe 

2 (9 + ADls ~ ^yPZ^yS + Dls) 



(116) 



f 4kTT 

=a|^ + 21og 



+ log 



(117) 



SV^VD^Vi + Dls 



(118) 



The caustic is a four-cusped astroid with the same an- 
gular extension along both axes, as can be explicitly seen 
transforming the above expressions to coordinates cen- 
tered on the caustic (see Ref. [l^l)- The outstanding 
feature that emerges from these expressions is that the 
drift of the caustic diverges logarithmically as the source 
approaches the horizon. As can be seen in Fig. [S^, the 
drift is always negative (clockwise as seen from the north- 
ern pole of the black hole) and grows linearly with the 
caustic order k. The drift tends to the asymptotic value 
(given by the first row of Eq. (|117p ) for large values of 
DlSj while it grows monothonically as Dls is decreased. 

The amplitude of the caustic does not present any di- 
vergences. As shown in Fig. [5)3, the amplitude increases 
linearly with the caustic order and tends to the asymp- 
totic value (given by the first row in Eq. (jllSp for large 
Dls- Dls is decreased up to the horizon it grows by 
a fixed amount, given by 

Rk{l)-Rk{^) = ^a^l^fil). (119) 

Taking Dls as a parameter ranging from 1 to -l-oo, 
we can trace the whole caustic hypersurface using Eqs. 
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FIG. 7: A 3-dimensional view (in Boyer-Lindquist coordi- 
nates) of the whole caustic tube for k — 3, a — 0.1 and 
equatorial observer no = 0. The sphere is the horizon of the 
black hole. 



FIG. 6: (a) Drift of the center of the caustic from the optical 
axis as a function of the source distance; from the bottom up, 
the curves are for fc = 2, 3, 4, respectively, (b) Semiamplitudes 
of the caustics with k — 2, 3, 4 (from bottom to top) as a 
function of the source distance. Both plots are done for a = 
0.1 and equatorial observer fio ~ 0. 



(jllSp and (|116p . The result is shown in Fig. [71 where 
the caustic appears as a tube with the transverse section 
having the shape of a four cusped astroid. At large dis- 
tances the caustic tube keeps its angular extension fixed 
and thus covers a larger and larger area. Close to the 
horizon, the caustic tube winds around the black hole 
indefinitely. A similar picture has already been done by 
Rauch & Blandford Q by numerical techniques for the 
caustic of order A: = 1. Our plot is entirely based on our 
analytical formulae (jllSp and (|116p . which are valid for 
arbitrary order except k = 1. Our study is thus comple- 
mentary to that of Ref. [3] . The origin of the logarithmic 
divergence in the caustic angular position can be traced 
back to the divergence of the integral I2 in Eq. (1751) 
(see also Eq. (|B4[) in Appendix |B| . Indeed, the inte- 
grand contains a factor A~^, which diverges linearly as 
the integration variable r approaches the horizon. Since 
Dls is the lower bound of that integral, I2 must diverge 
logarithmically as D^s tends to the horizon. The diver- 
gence of I2 has no effect in Schwarzschild, since it appears 
multiplied by the black hole spin, but as soon as a 7^ 0, 
the logarithmic divergence is transferred to the azimuthal 



shift 4>o — 4>Sj so that photons emitted by a source very 
close to the horizon must perform several loops before ex- 
iting. The latter argument is completely independent of 
our perturbative expansion in the black hole spin, prov- 
ing that the logarithmic divergence is not an artifact of 
our perturbative framework. Moreover, the divergence 
of I2 is not even a product of the strong deflection limit, 
since it is still there whatever the values of the constants 
of motion J and Q. Therefore, it seems plausible to us 
that even the primary caustic tube (with fc = 1), which 
is not included in our treatment, should show a similar 
behavior, winding an infinite number of times before en- 
tering the horizon. This seems not to be observed in the 
numerical analysis by Rauch & Blandford Q, where the 
primary caustic tube always appears to perform a finite 
number of loops before plunging into the horizon, except 
for the extremal case a = 0.5. It is anyhow difficult to 
believe that the divergence is compensated by any of the 
remaining terms in Eq. (|76p for intermediate values of a. 



Our formulae (|115p and (|116p represent the transverse 
section of the caustic tube at fixed -Dls- One might 
be interested to a different transverse section, e.g. the 
section at fixed (ps- This is particularly interesting to 
study the approach of the caustic tube to the horizon, 
when the drift is large and D^s is very close to 1. In this 
approximation, we can find the following expression for 
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any fixed value of 05 



MS = (-l)Vo ± Rk{DLs = 1) VI - ^o(l - ^')'^'' 

(120) 

16 



Dls =l + ye-^^{l 



2cfc - (2cfc + 7%/3) cos^ d i , (121) 



where A</) = 4>s — ^4>k{DLS = oo) is the difference be- 
tween the chosen value of (j)s and the asymptotic drift 
of the caustic center. This expression is strictly valid for 
large values of At/), but since the caustic drift sensibly 
increases only when Dls is very close to 1, it is sufficient 
that A(j) > 0.1 in order to validate this expansion. Fig. [H] 
shows the transverse sections of a caustic tube obtained 
at different values of A0. As the drift is increased, the 
caustic becomes thinner and thinner while it approaches 
the horizon. 



C. Gravitational lensing near caustics 

The higher order images of ordinary sources like stars 
or X-ray binaries are usually very faint, except for the 
event of a caustic crossing. Therefore, although in prin- 
ciple it is possible to analyze the lens equation in the 
general case, it is much more interesting to study the 
gravitational lensing of a source in the neighborhood of 
a caustic. This case is certainly the most relevant for ob- 
servations and is worth a complete and detailed analysis. 



1 . Position of the images 

The position of a source near a caustic can be expanded 
in the following way 



is = (1 - k)7r - A(l)k 



(122) 



(123) 



with dfis and S(j)s being of the same order of the caustic 
extension Rk, thus weighing as in the perturbative 
expansion. 

Correspondingly, the solutions of the lens equation will 
be very close to the critical points. We can thus use the 
following expansion for ip 



with Sil> being of order a^. 



(124) 



Using these expansions in Eq. (|105p . Si/j can be ob- 
tained as a function of the source position and ^ 



3V3^/d^{2Dls - 3)(3 + Dls?^^ 
{(3 - 2^ll - |2)(3 + DlsK^Dls - 3) 
9 - 2Dls - 2{2Dls - 3)i^D^ - ^i + DLs? 
47?Ls(3 + DLs){'iDLs - 3) (5 - 1') + lOS^Lsf } 



(125) 



whereas ^ is determined by the equation 



(126) 

in which S* = ±1 is inherited by the sign ambiguity of the 
(■0,^) parametrization (see [24I for more details about 
the resolution of the sign ambiguity). Note that the lens 
equation formally remains the same as in the D^g = 00 
case, though the position and the extension of the caustic 
have changed. The ^ equation (|126p can be easily put in 
the form of a fourth degree polynomial equation. It ad- 
mits two solutions if the source is outside the caustic and 
four solutions if the source is inside the caustic. The solu- 
tions so obtained satisfy the original equation (|126p with 
one definite choice of the sign S, which is then univocally 
determined for each image. 

Once we have the position of the image in the (V'jO 
space, it is straightforward to write the position of the 
image in the observer sky. To this purpose, it is impor- 
tant to note that ^ is known through Eq. (|126[) to zero 
order only. So, the position of the image on the observer 
sky is consistently determined to zero order as 



DoL0i= -^a^ + ek) 



(127) 



Dol02^ S^y/l~^il + ek), (128) 

with ek given by Eq. (|7T1) . The image appears on the 
critical ring of order fc at a position angle arcsin^ and 
half-sky determined by the sign S. 



2. Brightness of a lensed image 

In typical gravitational lensing studies, the change in 
the apparent brightness of the source is simply given by 
the geometrical magnification, defined as the ratio of the 
elementary angular area of the image and the angular 
area of the source as it would be seen without any grav- 
itational lensing. 
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When the source is far from the lens, this definition 
makes sense, since the background metric is asymptoti- 
cally flat and we can speak of a source without the lens 
by simply re-interpreting the source coordinates as coor- 
dinates in the asymptotic Minkowski metric. This proce- 
dure loses any meaning when the source is deeply within 
the gravitational field of the black hole. Moreover, the 
frequency of the emitted photon does not coincide with 
the frequency detected by the observer because of grav- 
itational redshift. Conservation of the photon number 
warrants that the quantity I /v'^ is conserved on a bundle 
of light rays, with I being the specific intensity, defined 
as the energy dE crossing a surface element dA pointing 
an angular area dO in the time interval dt and frequency 
interval dv. 

In order to build a simulated lensed image of a source 
close to a black hole, one just needs to find the position 
of the image for any point belonging to the source and 
determine the specific intensity observed at that point. 
The position of the higher order images for any given 
source position can be read from Eqs. (|127p and (|128p . 
whereas the specific intensity measured by the observer 
is related to the specific intensity emitted by the source 
through the relativistic relation 



I --^I 



(129) 



where the redshift factor can be calculated as usual as 

Po 



(130) 



with — g^ux'^ being the momentum of the photon 
and being the 4-velocity of the emitting particle. In 
stationary spherically symmetric and in Boyer-Lindquist 
coordinates for the Kerr metric, dt is a Killing vector and 
thus Po is a conserved quantity (we have put it to 1 by a 
choice of the affine parameter) . 



So, for any point of the source we can find the location 
of the corresponding high order images and if any model 
provides us the specific intensity of the source at that 
point, we can calculate the specific intensity as seen by 
the observer. 

In the case in which the source is transparent, it can 
be conveniently characterized by its volume emissivity 
ju{x^ ,p^), defined as the energy dE in the frequency 
interval dv^ emitted by the proper spacetime volume 
y/—gd'^x centered on in a solid angle dfl centered on 
the direction given by the vector p^. 

If we want the specific intensity measured by the ob- 
server in the sky direction (6*1, 62), we just have to trace 
back the null geodesic reaching the observer with such 
a direction and sum up the contributions given by all 
source elements along this geodesic. We thus have [35| 



Io = J '^Mx^,it.)dl 



prop-; 



(131) 



where dlprop is the geodesic line element as measured in 
the emitter frame. 



IV. CONCLUSIONS 

Year after year, our knowledge on the environment sur- 
rounding the supermassive black holes is growing at a 
higher and higher rate, thanks to the surprising develop- 
ment of the technology related to high resolution obser- 
vations. Angular resolutions of the order of the microarc- 
second are now reachable in the radio band and sooner 
or later will be reached in the sub-mm and X-ray band. 
At the same time, the interferometric observations in the 
infrared at the Keck telescopes and at VLT are unveil- 
ing a very rich stellar enviroment around the Galactic 
center, which is drawing more and more interest. High 
resolution observations are the premise for the discovery 
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of possible signatures of general relativity effects from su- 
permassive black holes, which would open a new era in 
the understanding of gravitational physics. 

Within this context, the comprehension of the propa- 
gation of photons in a strong field environment is of cap- 
ital importance. Many numerical codes partially exploit- 
ing the analytical solutions of Kerr geodesies are available 
and have been used to build simulated images of the black 
holes. At the same time, analytical studies have been 
complenientarily developed to conquer precious insight 
about the mathematical structure of the lens mapping in 
this very special gravitational lensing framework. 

The strong deflection limit allows to study gravita- 
tional lensing in the extreme situation of photons trav- 
elling very close to the unstable circular orbit around 
the black hole. Such photons give rise to an infinite 
sequence of additional images which contribute to the 
total flux received by the observer by a non-negligible 
amount (see Ref. 3 for an estimate of their relevance) . 
Therefore, the study of such images is far from being 
a mere academic exercise, but acquires a striking impor- 
tance by the fact that these images carry invaluable infor- 
mation about the strong gravitational fields around the 
black hole. Whereas complete numerical studies of these 
images are difficult because of the extreme accuracy re- 
quired to follow photons travelling around the unstable 
circular orbit, analytical studies benefit from the great 
simplification introduced by the strong deflection limit. 
Therefore, higher order images represent a unique win- 
dow where we can confront simple analytical results from 
General Relativity with observations. 

In this paper we have removed the traditional limi- 
tation of the strong deflection limit to sources very far 
from the black hole. We have thus analytically explored 
extreme gravitational lensing of sources close to a black 
hole for the first time. There is no limitation to the va- 
lidity of our results, which are applicable even to sources 
just outside the horizon. 

In the spherically symmetric case, we have shown that 
the same formulae for the deflection of the photon can be 
applied both to sources outside the photon sphere and to 
sources inside the photon sphere, whose images appear 
inside the so-called shadow border. We have specified 
our formulae to the Schwarzschild case in order to test 
the validity of the strong deflection limit throughout the 
range of source distances. 

For the Kerr black hole, the only modification comes 
in the resolution of the radial integrals in the geodesies 
equations. This brings to a modification of the size and 
shape of the critical curves. We have obtained a com- 
plete analytical description of the caustic tube from infi- 
nite distances up to the horizon, showing that the caustic 
tube winds indefinitely around the black hole because of 
a logarithmic divergence in the azimuthal geodesic equa- 
tion. We have also updated the gravitational lensing of 
sources near a caustic. 

The formulae contained in this paper can be applied to 
physically motivated models of sources around a super- 



massive black hole, where they can be used to calculate 
the shape and the brightness of the higher order images. 
We leave this interesting task to future work, contenting 
ourselves with the complete analytical derivation of all 
relevant formulae here. 
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APPENDIX A: DIFFERENT PERTURBATIVE 
PARAMETERS FOR THE STRONG 
DEFLECTION LIMIT 

Recently Iyer and Petters [3l| have rewritten Darwin's 
formula for the defiection angle of a Schwarzschild black 
hole in the strong defiection limit in terms of a new pa- 
rameter b' , defined as 

6' = 1-^. (Al) 

u 

Recalling the definition of Darwin's perturbative param- 
eter e (|19p . we get the simple relation 



Of course, if we stop at the lowest order in the expan- 
sion, the defiection angle can be indifferently expressed 
in the equivalent forms 

a = - loge + log[216(2 - Vif] - tt + 0(e) (A3) 
= - log 6' + log[216(2 - V3 f] - vr + 0{b').{AA) 

Both formulae are correct to lowest order, the differences 
being stored in the higher order terms in the perturbative 
expansions. Iyer and Petters have found that the higher 
order discrepancy between these two formulae and the 
exact deflection angle is significantly smaller for the b' 
formula than for Darwin's one [3l|. Of course, by spec- 
ifying an invertible function f\{€) such that /a(0) — 
and f'x{0) 7^ 0, one can always define A = f\{e) as a new 
perturbative parameter. By a suitable choice of f\ one 
can make the higher order corrections vanish up to an 
arbitrary order n. However, this normally corresponds 
to a more and more complex form of f\, which spoils the 
advantages of the perturbative expansion. The choice of 
the perturbative parameter is arbitrary within the men- 
tioned constraints on />, but is normally driven by some 
physical quantities that can be easily expressed in terms 
of the perturbative parameter. 

In the case of gravitational lensing in the strong de- 
fiection limit, this point can be made clearer once we 
construct the formulae for the position of the images in 
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the two perturbative frameworks. If we use the b' for- 
mula in the lens equation, we trivially obtain b'^ — e„ for 
the n-th image, with e„ always given by the expression 
([5^ . The difference between the two frameworks actu- 
ally emerges from the expression of the angular position 
in the observer sky 6 in terms of the new parameter, 
which reads 



hole spin a. Actually, the integral I2 is already multiplied 
by a in Eq. l[76|) . so that it is sufficient to stop at the 
first order in its expansion. The results are 

/i = -ailogS + bi+ci{DLs)+ci{DoL) (Bl) 

/2= -a2l0gS + b2+C2{DLs)+C2{DoL), (B2) 



OL 



+ = 



1 - b' 



(A5) 



If we plug 6^ = e„ directly into Eq. ()A5|) we have a new 
formula for the position of the images 



On 



l-er. 



(A6) 



which can be compared to the classical formula by Dar- 



(A7) 



In both formulae e„ is given by Eq. ()54p and is just a 
function of the source position, with no memory of the 
perturbative framework used. 

If e„ is small, the two formulae are very close each 
other, the difference being order and thus completely 
negligible. If an accuracy is necessary, both formulae 
must be complemented by their respective higher order 
terms. Then the difference will be in the third order and 
so on. 

On the other hand, in intermediate situations, in which 
e„ is not small, Eq. (jA6p does significantly better than 
Eq. (jA7P and can be used to extend the range of valid- 
ity of the first order expansion of the Strong Deflection 
Limit, without resorting to the second order terms. 



APPENDIX B: ADDITIONAL TERMS IN THE 
RADIAL INTEGRALS 



As anticipated in Section UlI A[ the dependence on the 
source and observer distances entirely comes from the 
resolution of the radial integrals ([M]) and ((99)) . These 
integrals can be performed using the general tools of Sec- 
tion In] and then expanded to second order in the black 



with the coefficients ai, 61, a2, and 62 being unchanged 
with respect to the appendix of Ref . [22] ■ To the second 
order in a, the function ci reads 



ciW = y log 



(2-f Va) 



3r - yS+T 



3r + V3+T 



3\/3V3T7(2r - 3) 
W [27(2r - 3)2^^(3 + rf/^ ^ 
■ |2(2r - 3)(3 + r) [2(2r - 3)(Vf - Vf+V'fil - f ) 
+ I4:r+I0r^^ -9{l-i^)+8rfil -216re^|. (B3) 

To the first order in a, the function C2 is 



C2(r) 



02 
2 

log 
2a 



log 



/3f- vsT^ 



/3f -I- 



2V^+ x/3T7 



3(2V^- V3T7) 



+ [V^V3T7(2r-3)]" 
• [2(2r - 3)(Vr - VsTrf + 18r - 9] . (B4) 

It can be easily checked that both functions vanish as 
their arguments go to infinity. Moreover, if r — > 3/2, the 
divergence in a is compensated by the vanishing of S. As 
r < 3/2, it is also easy to see that we must have (5^ < 0, 
analogously to the spherically symmetric case. The sec- 
ond logarithm in C2 diverges at r = 1. The implications 
of this divergence are discussed in Section IIII Bl 
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